Universal rate constants for reactive collisions of ultracold molecules 
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A simple quantum defect model gives analytic expressions for the complex scattering length and 
threshold collision rates of ultracold molecules. If the probability of reaction in the short-range part 
of the collision is high, the model gives universal rate constants for s- and p-wave collisions that are 
independent of short-range dynamics. This model explains the magnitudes of the recently measured 
rate constants for collisions of two ultracold ^'^K^^Rb molecules, or an ultracold '"'K atom with the 
'"'K'^^Rb molecule [S. Ospelkaus et al, Science 327, 853 (2010)]. 
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Research utilizing ultracold atoms has become estab- 
lished as a major forefront multidisciplinary area involv- 
ing diverse areas such as atomic and molecular physics, 
quantum optics, condensed matter physics, and quan- 
tum information. Extending such research with ultra- 
cold molecules will open up a number of new opportuni- 
ties [H . Development of sources of ultracold molecules 
is progressing rapidly Q . It is essential to understand the 
elastic, inelastic, and reactive collisions of such molecules 
in order to control and utilize them effectively. A re- 
cent experiment reported the first evidence of ultracold 
"chemistry" and measured the reactive collision rates of 
fermionic ^°K*^Rb molecules between 250 and 900 nK Q- 
Since molecules are much more complex than atoms, it is 
necessary to develop models of threshold molecular col- 
lisions that are adequate to account for their complexity 
yet simple enough for understanding experimental data. 

Here wc propose a simple, yet general, model for 
threshold molecular collision rates based on the separa- 
tion of the effects of the long- and short-range parts of 
the intermolecular potential for a given scattering chan- 
nel 0. Our fully quantum mechanical model uses the 
analytic framework of generalized multichannel quantum 
defect theory (MQDTMH] based on the specific version 
pioneered by Mies [SMlOj. The reactive or inelastic part 
of the short-range collision is characterized by a dimen- 
sionless "quantum defect" parameter < y < 1, which is 
related to the probability of irreversible loss of incoming 
scattering flux from the entrance channel due to dynam- 
ics at short range. The role of the long-range potential 
is to determine how much of the entrance channel wave 
with very large de Broglie wave length is transmitted 
to short range to experience such loss dynamics. Our 
MQDT model generalizes the model of Ref. [HI, which 
has been successfully used to understand the magnitude 
of inelastic collision rates involving ultracold atoms. 

We find universal rate constants for s- and p-wave col- 
lisions depending only on the long-range potential when 
y ^ I, which corresponds to unit probability of reaction 
at short range so that no flux is reflected at short range 



back into the entrance channel. In contrast, if y ^ 1, the 
rate constants depend strongly on the scattering length 
of the entrance channel and are not universal. 

We consider collisions of two particles interacting via 
the van der Waals (vdW) potential at long-range. These 
can be S'-state atoms, or molecules in the rovibrational 
and electronic ground state. The scattering channels 
|a) = \aia2)\irn) are defined in terms of their respec- 
tive internal states oi and 02 , and the partial wave quan- 
tum numbers £m. In general, the elastic rate constant 
JC°^ and the rate constant lO^ for total inelastic or reac- 
tive scattering are given by the diagonal elements Saa of 
the S'-matrix for channel a. Expressing Saa = e^"^° in 
terms of an energy-dependent complex scattering length 
a-aik) = aa(fc) — iPa{k) [l2j, defined analogously to the 
s-wave energy-dependent scattering length [ISj, , 

~ ^ t&\nja(ka) 11 Saaikot) ^-y^ 

ka ika 1 + Saa{ka) 

gives the rate constant contributions from channel a, 

^t{E) = ga^ |1 - SaaiEf = 2ga — \~aa{ka)?fM , 



l^aiE) = ga^ (1 - \Saa{E)\^) = -/3„ (fc„ )/„ 



(2) 
(3) 



where fc^ = 2^(£' — Ea)/h'^ with E denoting the total 
energy, fj, the reduced mass, and Ea the threshold energy 
of the channel a. The factor ga = ^ except that ga = 
2 when both particles are identical species in identical 
internal states, ai = 02; £ is restricted to being even 
(odd) in the case of identical bosons (fermions). The 
function 



fa{ka) — 



1 



1 + kl\ha{kaW + 2kaPaika)' 



(4) 



has the property that < fa{ka) < 1 and fa{ka) 1 
as ka 0. The latter is true when both conditions 
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ka\aa{ka)\ ^ 1 and kapa{ka) 'C 1 are met. We stress 
that Eqs. (HI)-© are exact [12|, and they apply either to 
coupled channel (unitary) and complex potential (non- 
unitary) models. 

We introduce a single channel model with a complex 
potential to represent the long range potential and short 
range loss dynamics for channel a, where for simplicity 
of notation we drop the implied channel index a: 



Uiir) = Vir) + 



h^£{£+l) . 7(r) 



(5) 



and V{r) has the vdW form at long-range: V{r) = 
—Cq/t^ for r > Rq, where Rq denotes the range of 
short-range forces (e.g. the exchange interaction). The 
vdW potential is characterized by the length Rq = 
^ (2fj,CQ/h'^Y^ or the closely related length a = 
47ri?6/r(|;)^ [la]- The imaginary part 7(r). which is as- 
sumed to vanish beyond Rq, simulates all short-range 
coupling at r < i?o to exoergic non-threshold exit chan- 
nels that result in loss from the entrance channel. We 
assume that the kinetic energy release in loss channels 
is much larger than any exit barriers. The short range 
V{r) for r < Rq should be viewed as a pseudopoten- 
tial which determines the phase of the wave function for 
r > Rq. This phase determines the "background" s-wave 
scattering length a for the potential V{r), which can in 
general be different for different £m in the case of general 
anisotropic potentials. 

Decomposing the wave function into real and imagi- 
nary parts, \E'(r) = F(r) + iG{r), gives the two channel 
Schrodinger equation 



5^ 

Q^2 



> 



W(r)) *(r) = 0, 



(6) 



for the wave function $(r) = {F{r),G{r)} with a non- 
Hermitian interaction matrix 



W(r) 



2 




(7) 



We solve the two-channel problem in the framework of 
MQDT 0-01 , which separates the effects of the long- and 
short-range parts of the potential. We adopt the nota- 
tion of Mies 0, and introduce a quantum defect ma- 
trix Y(£:) = {{0,-y{E)},{y{E),0)}}. The antisymmet- 
ric form of Y results from the non-Hermitian potential 
W(r). 

The wave function for r > Rq is 



*(r) = /(r)I +g{r)Y{E) 



(8) 



Here, A is some constant vector, I is the identity matrix, 
and /(r), g(r) are solutions of the Schrodinger equation 
in the real part of the potential Vi{r) that have "local" 



WKB-likc normalization at short distances (s. 16| 

f{r,E) 9i k{r)-^/^ sin I3{r), 



){r,E) = k{ry^/'^ cos I3{r), 



r > Ro, (9) 



where k{r) = ^J2n{E — 'Re\U i{r)]) / h is the local wave 
vector, and /3(r) = ^^Axk{x) is the WKB phase. Typ- 
ically, Ro Re, and for s- or p-wave ultracold collision 
one can assume that y{E) = y is independent of E. The 
two parameters a and y completely represent the effects 
of the complex short-range dynamics on the wave func- 
tion in channel a for r > Rq. 

Some intuition about y is gained from writing *I'(r) at 
small distances, Rq < r ^ i?g, using the WKB-like form 
in Eq. © for /(r) and g{r): 

exp [—i J^k{x)dx] f ^ ^ y\ ^-'^P [* J^k{x)dx] 
l + y 



(10) 

The first term represents the flux of incident particles, 
whereas the second term gives the flux reflected from the 
short-range potential. Hence, for y = 1 there is no outgo- 
ing flux, whereas for y not unity there is some back reflec- 
tion. Finally, for y = (no losses), the incident and re- 
flected fluxes arc equal, and we recover the standard scat- 
tering wave function ^(r) ^ k~^^'^{r)sin^J^k{x)dx'\ = 
f{r). In the following we will assume < y < 1 

Applying the standard MQDT formulas relating the 
quantum-defect matrix Y to the scattering S-matrix , 
we calculate a{E) assuming kR^ <C 1 



tanC(-B) 



(11) 



oo, (12) 



i -I- y tanA(i?)^ 

Here, C{E) and tanA(£;) are the MQDT functions [3 
that connect the solutions with short- and long-range 
normalization: /(r) = C-^{E)f{r), g{r) = C{E)[g{r) + 
tan A(iJ)/,;(r)], where the scattering solutions in the 
asymptotic region are defined as follows: 

/(r, E) = sin {kr - ^7r/2 + C) /Vk, 
g{r, E) = cos {kr - £t:/2 + ^ /Vl, 

where ^(fc) is the phase shift. The functions C{E), 
tan ME) and (,{E) are found from the analytic vdW the- 
ory [l8|. The small- fc behavior for an s wave is 

C-2(i;,£ = 0)^fca(l + (.s-l)2), (13) 

tanA(i;,£ = 0) ^ 1 - s, (14) 

and tan^(£',^ = 0) ^— ^ — fca, where s = a/a is the di- 
mensionless scattering length. For a p wave: 

^2l + (s-l)' 



C-'^{E,£ = 1)'^ 2kai{kay 

tan\{E,i^ 1) s/(s- 2), 
tan£_{E,i = 1) 2kai{kaf 



{S-2Y 



(15) 

(16) 
(17) 



where cii = ar(i)V(1447r2r(|)2) w 1.064a. 

For small y Eqs. ([3]) and pT|) predict that the loss rate 
scales as /?(£') « yC~'^{E). This can be interpreted 
as an energy-insensitive probability y of the short-range 
reaction multiplied by the fmiction C~^{E), which pro- 
vides the proper scaling of the wave function amplitude 
between the short-range and asymptotic zones. 

A key result of this paper is Eq. (fTTj) for d{E) for ar- 
bitrary y. Equations (fT3 |) - (fT7| yield the small-fc analytic 
limits of Eq. (fTTj). valid when k\d\ <C 1 and kR^ ^ 1: 



di^oik) = a + ay- 
de=i{k) 



i + y{l- s) 
2ai(fca) 



2 2/ + i(s 



1) 



ys ■ 



2) 



(18) 
(19) 



These expressions parameterize the elastic and inelastic 
collision rates in terms of the two dimensionless param- 
eters s and y, namely, an entrance channel phase and a 
short-range interchannel coupling strength, plus the van 
der Waals parameter a. In the special case of unit prob- 
ability of short range loss from the entrance channel, 
?/ — !■ 1, we find the universal result that a is indepen- 
dent of s and y, and depends on d only. Specifically, 
Q.£=o — Pi^Q = d for an s-wave and ai=i = — ai(fca)^ 
and f3e=i — > ai(fca)^ for a p-wave. This gives the uni- 
versal rate constants determined solely by the quantum 
transmission of the long-range potential: 



h 



(20) 



ICgli = ^c7g-kd\{kdY = 2ag-di{kdf (21 



valid for fci?g ^ 1. The factor cr = 3 for ^ = 1 for 
a rotationless molecule in the universal limit, since all 
three m components have the same Cq and contribute 
equally. In general, a sum of the different contributions 
from each m would need to be taken. 

We have verified Eqs. (PU)) and (|2ip by numerical cal- 
culations with a complex potential. These equations also 
reproduce the numerical model of Orzel et al [llj for s- 
and p- waves, for which y = 1 is an excellent approx- 
imation for ionizing collisions of Xe metastable atoms. 
Ref. 
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found the ionization rate constant for s-wave 
collisions to be independent of isotopic mass, that is, in- 
dependent of s, which is expected to be quite different 
for different isotopes The value of ICf^Q predicted 



by Eq. (f20| agrees with the measured value of Ref. [11 1 
within experimental uncertainty. Hudson et al. [l9| also 
found that a numerical implementation of the model of 
Ref. [m explained their vibrational relaxation data for 
RbCs molecules, although the error bars were quite large. 

In general, y is not unity, and there can be reflected flux 
from short-range back into the entrance channel, so that 
d is sensitive to both s and y. Figure [T] shows the complex 
scattering length d plotted in the {a, 13} plane as y and 
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FIG. 1: Real and imaginary parts of af=o/a (left panel) and 
(a£=i/a)(A:a)~^ (right panel) as — ^ for different values of 
the loss parameter y (set of circles) as the phase parameter 
s = a/a, color coded according to the value of arctan(s)/7r, 
varies over its full range of —oo < a < +00. The blue-red 
boundary occurs where Ue{r) has a bound state at .E = 0. 



s vary for s- and p- waves. The Figure shows d lies on a 
circle for a fixed values of y. For s-waves these circles are 
centered at {a, d{^ + ^)}, and have radius a(l — y^)/(2y), 

while for p- waves they are centered at {— ai,ai(| -I- ^)} 
and have radius ai(l — y'^)/{2y). The points at the top 
(bottom) of the circles correspond to the maximal (mini- 
mal) loss: l3i=o = d/y and (3e=i = di/{yk'^d'^) {/3e=o = yd 
and = yai/(fc^a^)), and they are realized at s = 00 
and s = 2 [20| (s = 1 and s = 0), for the s- and p-wave, 
respectively. Hence, counterintuitively, the largest loss 
rates (/3 ^ d) occur only for small values of the reac- 
tion probability y, but only when s > 1 is large enough 
that there is a near-threshold bound state that allows 
the C(-B)^^ function to build up a large amplitude of 
the short-range wave function. 

Our theory explains recent experimental data on the 
reaction rate coefficients of a gas of ultracold ''^^K^^'Rh 
molecules in their vibrational and rotational ground 
state These molecules are fermions with 36 distinct 
states of nuclear spin. The experiment could prepare 
the molecules in either the same or different spin levels, 
thereby allowing measurement of both s- and p-wave re- 
active collision rates. Fig. ^ compares the experimental 
data with two different predictions of our MQDT model. 
The dashed line shows the predictions of the universal 
model at low fc, where ^^Lol-^) ^ 4(/i//i)a is indepen- 
dent of temperature T and 



/Cti(T) = 



r(i/4f 
r(3/4)2 



f^— = 1513a^-?— 



(22) 



varies linearly with T. With a = 118(3) ag from two 
ab initio calculations of Cq [21|, , assuming a 5 % un- 
certainty, ICf^i(T)/T = 0.8(1) X 10" 



pared to the measured value of 1-1(3) x 10^^ cm'^s^^K^^ 
for the lowest spin state 

While the universal p-wavc ICf^^ agrees with the data 
within mutual uncertainties, the s wave /C'/__n — 0.8 x 
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FIG. 2: Loss rate constant K.^^ of KRb molecules versus tem- 
perature. Experimental data for a single component gas 
of molecules in spin states \F — —4, Mp = ^) (blue squares), 
\F = —4,Mf = |) (red triangles), and for 50/50 mixture of 
these two spin states (green dots), are compared with predic- 
tion of MQDT in the universal regime {y — 1, dashed lines), 
and in non universal regime {y < 1, solid lines). In the latter 
case the values of y, and s parameters for i = and i = 1 are 
determined by fitting to the experimental data. 



10^ cni^/s is a factor of 2 smaller than the measured 
1.9(4) X 10"^" cm^/s. This departure from universal- 
ity may yield information about the short range dynam- 
ics. We have fitted the formula ([3|) with 5 and (3 from 
Eqs. (fT8|) and (fT9|) to both the s- and p-wave experi- 
mental data to estimate values of the y and s parame- 
ters, assuming two distinct s parameters for £ = and 
£ ~ 1. Since the equations for the complex scattering 
length are quadratic in s, this yields two sets of solutions: 
(i) y = 0.397, Si=o = 10.3, Si=i = 3.43; (ii) y = 0.420, 
Si=o = -48, si=i = 3.43, The sohd lines in Fig. ^ show 
that agreement can be obtained with the data in this 
way, implying that a significant fraction of the incoming 
flux may be reflected back into the entrance channel. 

Our model predictions can be tested by other exper- 
iments. For example, if y is not unity, unlike the case 



studied by Ref. [Ill ] , there would be some dependence on 
isotopic mass, since s will vary with mass [l5l.[l6j. Thus, 
other isotopic species would depart from the universal 
rate constants differently than the "^^K^^Rb isotope. In 
addition, experiments with an electric field turned on to 
induce a molecular dipole should have short range reac- 
tions governed by the same s and y parameters, if the 
asymptotic solutions of the dipole problem are matched 
to the solutions in the vdW zone Rq < r < Rq. Thus, 
collision rate measurements in an electric field could test 
universality and the role of short range dynamics. 

Equation (|21|) predicts a universal rate constant of 
1.1(1) X 10~^° cm^/s for reactive s-wave collision of 
atoms with ''°K*^Rb molecules, using the Cg constant 
from Kotochigova [2lj. This is quite close to the mea- 
sured value, 1.7(3) x 10~^° cm^/s, of Ospelkaus et al [1]. 

Our MQDT parameterization provides a simple frame- 



work for characterizing and understanding threshold 
molecular collisions. The MQDT method could readily 
be generalized to include the effects of molecular reso- 
nance states, threshold exit channels, finite electric fields, 
or reduced dimensional geometry. Collisional resonances 
should be suppressed when j/ — > 1 but should be promi- 
nent if 2/ ^ 1. A MQDT framework may also offer a way 
to incorporate quantum dynamics into simpler general- 
ized quantum threshold Langevin models of cold molec- 
ular collisions in electric fields (23| . 
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